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ABSTRACT 

A flow diagnostic technique based on the molecular 
Rayleigh scattering of laser light is used to obtain 
dynamic density and velocity data in a high speed flow. 
The technique is based on analyzing the Rayleigh 
scattered light with a Fabry- Perot interferometer used in 
the static, imaging mode. An analysis is presented that 
established a lower bound for measurement uncertainty of 
about 20 m/sec for individual velocity measurements 
obtained in a 100 psec time interval. Software and 
hardware interfaces were developed to allow computer 
control of all aspects of the experiment and data 
acquisition. The signals from three photomultiplier tubes 
were simultaneously recorded using photon counting at a 
10 kHz sampling rate and 10 second recording periods. 
Density and velocity data, including distribution functions 
and power spectra, taken in a Mach 0.8 free jet are 
presented. 

INTRODUCTION 

In this paper we describe a flow diagnostic based 
on the molecular Rayleigh scattering of laser light for 
use in NASA test facilities. This diagnostic will provide 
dynamic aerothermodynamic data that are not presently 
available. These data will be important in studies such 
as the investigation of growth and decay of turbulent 
fluctuations. The effort is part of the non-intrusive 
instrumentation development program supporting 
propulsion research at the NASA Glenn Research 
Center. This work is directed to the measurement of 
fluctuations in flow velocity, density, and temperature 
for jet noise studies. One of the main objectives in jet 
noise research is to identify noise sources in the jet and 
to determine their contribution to noise generation. In 


particular, researchers have focused on the correlation 
of fluctuations in flow parameters with far field noise 1 . 

A variety of nonintrusive, laser based flow 
diagnostics (Rayleigh scattering, Laser Doppler 
Velocimetry (LDV), Particle Imaging Velocimetry (PIV), 
Laser Induced Fluorescence (LIF)) are routinely being 
applied for time average and instantaneous planar 
measurements of velocity, density, temperature, and 
species concentrations. However, an important class of 
measurements, high frequency response dynamic 
measurements of flow parameters, is not addressed by 
current laser diagnostics. This type of time history data is 
needed to determine, for example, density and velocity 
spectra, density-velocity correlations, and two-point 
correlations. Although LDV can typically achieve data 
rates of a few tens of kHz, it is difficult to achieve much 
higher rates without introducing larger amounts of seed 
material into the flow. In practice, LDV normally 
provides mean velocity and a measure of turbulence 
intensity. Furthermore, the measurements are random 
in time, which makes it difficult to obtain time history 
data needed for power spectra. In turbulent flows, LDV 
measurements are beset by a variety of so-called biasing 
errors, caused by correlations between the measurement 
rate and flow properties. Planar techniques, such as PIV 
and planar Rayleigh scattering, provide a large number of 
simultaneous measurements in a plane, but are generally 
limited to low sampling rates, determined by the pulse 
repetition rate of the laser and by the time needed to 
transfer image data from the camera. New techniques 
are needed to provide nonintrusive, dynamic 
measurements that can provide data similar to that 
provided by hot wire anemometers. Since it is unrealistic 
to expect (at least at the present time) to make 
measurements at a large number of locations and at a high 
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sampling rate, it seems prudent to develop a laser 
diagnostic capable of point measurements at high 
sampling rates. 

Because the Rayleigh scattering technique under 
study is based on molecular scattering rather than particle 
scattering, no seed material need be injected into the flow. 
One difficulty with LDV and PIV is that they require the 
flow to be " seeded" with micron size seed particles to 
provide a sufficient concentration of scattering centers. 
The reliance on seed particles presents a number of 
difficulties, such as difficulty injecting a uniform cloud of 
seed, and providing seed material that can withstand the 
flow environment. For high temperatures, refractory seed 
materials are necessary. In addition, seed material can 
contaminate the facility, coating surfaces and windows. 
Facility engineers are often wary of the introduction of 
this foreign, often abrasive material, into their equipment. 
A further limitation of particle scattering methods is that 
the particles, although small, may not be able to follow 
large flow accelerations, thus introducing inaccuracy in 
the prediction of the gas velocity. 

These problems associated with particle scattering 
measurements are eliminated if molecular scattering is 
used, since the gas molecules that constitute the flow 
under study are used as the scattering centers. The 
simplest molecular scattering based diagnostic is 
Rayleigh scattering. The frequency spectrum of Rayleigh 
scattering is closely related to the velocity distribution of 
the scattering gas. The spectrum may be analyzed to 
determine temperature, density, and velocity. Density is 
simply proportional to the total scattered light; 
temperature is related to the width of the Rayleigh 
spectrum; and one component of velocity is proportional 
to the shift of the spectral peak from the frequency of the 
incident light. Because the spectral width is also a 
function of the molecular weight of the gases in the flow, 
knowledge of the gas composition is generally required; 
however, this is not a concern in the proposed work, 
which is directed toward air flows where the composition 
is well defined. In any case, velocity measurements, 
because they are determined from the frequency of the 
peak of the spectrum, are independent of the gas 
composition. Rayleigh scattering is particularly suitable 
for measurement of supersonic and hypersonic velocity 
where the mean molecular velocity (flow velocity) is 
larger than the random molecular velocity (temperature). 

Because of the relative simplicity of Rayleigh 
scattering based gas density measurements, they have 
been more widely used than velocity and temperature 
measurements. In our work, for example, time and phase 
averaged density measurements have been made in an 
underexpanded supersonic free jet in support of jet noise 
studies^. We have also demonstrated simultaneous 
multiple point density measurements for determining gas 


density fluctuation spectra, cross-spectra, and cross 
correlation functions in a low speed heated jet 3 . The 
more difficult problem of dynamic velocity and 
temperature measurements was addressed in two previous 
papers. The first paper 4 used direct imaging of the 
Rayleigh scattered light through a Fabry-Perot 
interferometer. One photomultiplier tube (PMT) was 
used to detect the total Rayleigh scattered light for density 
measurements (as is done in the present paper). Two 
additional PMT’s were used to determine the velocity and 
temperature. A disadvantage of this method, however, is 
that it required the detection system, including the 
sensitive Fabry-Perot interferometer, to be located close 
to the experiment. A second paper 5 described a system 
that used an optical fiber to transmit scattered light from 
the experiment to the Fabry-Perot Interferometer (FPI), 
which could be located at some distance from the 
experiment. This allowed the interferometer to be located 
in a controlled environment, where the temperature can 
be held constant and where acoustic noise and vibration 
can be minimized. This system also incorporated an 8x8 
multiple anode PMT that was used to detect the 
interference pattern formed by the FPI. However, the 
quantum efficiency of the 8x8 multiple anode 
photomultiplier was only about 5 %. This is low 
compared with the quantum efficiency of about 25 % that 
can be achieved with single PMT’s. In addition, both 
these systems used analog data acquisition equipment, 
which resulted in relatively high noise levels in the 
detected signals. 

In this paper we present an improved version of 
these systems. The present system uses three 25% 
quantum efficiency PMT’s to improve the detection 
efficiency. A new photon counting data acquisition 
system is used, which improves the signal to noise ratio of 
the detected signals by reducing the amount of noise due 
to electrical interference. Also, new computer software 
and hardware have been incorporated that greatly 
simplify the system operation. Some basic background 
material related to molecular Rayleigh scattering based 
diagnostic systems is first presented. The lower bounds 
for measurement of density, temperature, and velocity are 
given, assuming errors result only from Poisson noise 
generated in the photo detection process. This establishes 
the best possible measurements that could only be 
approached by use of an ideal instrument. We then 
describe the scheme used to measure dynamic velocity 
based on a Fabry-Perot interferometer and three PMTs. 
The lower bounds for velocity measurement uncertainties 
are calculated for this particular system. 

An experiment conducted to evaluate the technique 
on a subsonic flow (up to about 300 m/sec) is then 
described. Results are given for data sampled at a 10 kHz 
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rate. Examples of velocity time history and power 
spectrum of velocity fluctuations are given. 


THEORY 

Rayleigh scattering 

The spectrum for Rayleigh scattering from a low 
density gas has a Gaussian profile given by 


«/-/.)# =^-expi- 1 
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where f Q is the laser frequency and u is the mean gas 
velocity. The interaction wave vector is K = k-ko (with 
ko and being the wave vectors of the incident and 
scattered light), and a = (2 k77w) 1/2 is the most probable 
molecular speed (with k being Boltzmann's constant, m 
the molecular mass, and T the gas temperature). Note 
that the spectral peak is shifted by a frequency 
proportional to the component of the bulk velocity in 
the K direction. The spectral width is proportional to 
the square root of the gas temperature. It is convenient 
to introduce the velocity component u K = K* u IK, 
which represents the measured velocity component. 

The assumption of a Gaussian shaped Rayleigh 
scattering spectrum is only valid if 



where p is the gas pressure and rj is the shear viscosity. 
Collective effects of the molecular motions become 
important for higher density gases (y ~ 1) and a more 
detailed kinetic theory model, such as the Tend S6 
model 6 , is required to describe the Rayleigh scattering 
spectrum. Although the Tend spectrum does differ 
significantly from a Gaussian spectrum for the flows we 
are studying, similar results are obtained for the 
uncertainty analysis we conducted. We, therefore, used 
the Gaussian spectral model to reduce the 

computational time. 

Lower bounds for ideal measurements 

Since Rayleigh scattering is a relatively weak 
process, the uncertainty in the measurements often is set 
by the photon statistical noise (shot noise), which 
determines the lower bound on measurement 

uncertainty. For example, the variance in the number of 
photoelectron counts for a Poisson process is equal to 
the mean number of counts. Thus the lower bound for 
the relative uncertainty in the measurement of gas 
density p, is equal to the square root of the variance 
divided by the mean counts N R . The lower bounds for 
temperature and velocity uncertainties have also been 
evaluated for a low-density, one -component gas. The 


relative uncertainties caused by photon statistics for this 
case can be written 7 ’ 8 . 
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These relations provide a lower bound for the 
measurement uncertainties, which can only be attained 
if the shot noise is the dominant noise and if an ideal 
instrument is used. 

Consider an experiment with the following 
parameters: a laser (532 nm) with output power P 0 = 
5 W, air at STP (T = 293 K, p = 1 atm), probe volume 
length L x = 1 mm, f/4 collecting optics (i.e., solid 
collection angle Q - 0.05 sr), and efficiency factor £ = 
5 %. The rate of detected photons, given by 


jlj eP 0 nL x Xa(do\ . 2 
= 50 million counts / sec 


(3) 


In this equation, n is the gas number density, da/d£2 is 
the differential scattering cross section (6.13x10 " m7sr 
for air at 532 nm), x is the angle between the electric 
field vector of the (linearly polarized) incident light and 
the direction of the scattered light, h is Planck's 
constant, and c is the velocity of light. If we wish to 
obtain independent measurements at a 10 kHz rate, the 
total number of detectable photons in each period would 
be 5000. Equation 2 gives the lower bound for 
measurement of density, temperature and velocity for 
each period: 

a ^ = 1.4 %, = 2.0 % , a(u K ) = 4.1 m/sec 

p T y K) 

This shows that high sampling rate Rayleigh scattering 
measurements are at least feasible. It must be 
emphasized, however, that these values represent the 
best possible measurements given an ideal instrument. 
In practice, we are limited to instruments, such as the 
Fabry-Perot interferometer used here, which result in 
significantly higher uncertainties, as described below. 

Lower bounds for practical instrument 

Estimates of the measurement uncertainty for the 
technique described here, where the Rayleigh scattered 
light is analyzed with a planar mirror Fabry-Perot 
interferometer, are obtained by numerically calculating 
the Cramer-Rao lower bound 9 . The variance of the 
estimate of a parameter a, (e.g., temperature or 
velocity) is given by 
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V(a,.) = [r-'] (/ (4) 

where no summation over repeated indices is implied. 
For Poisson statistics, T is the Fisher information matrix 
with elements 

r _ y 1 (5) 

* . (n d ,) <**.■ da i 

where <N Dq >is the expected number of counts from the 
q lh photodetector. The summation is over the number of 
photodetectors, and the quantities <N Dq >depend on the 
details of the optical system and the flow parameters, as 
presented below in the Uncertainty analysis section. 

Optical configuration 

In order to perform an uncertainty analysis, we 
first describe the setup (fig. 1 ) used for the experimental 
work. Light from a 5W, 532 run, single-frequency, 
Nd: Vanadate CW laser was focused by lens LI 
(350 mm focal length) to a 150 pm diameter beam at 
the probe volume. The laser beam was terminated in 
light trap LT. Rayleigh scattered light is collected at a 
90° scattering angle and focused by L2 (two 160 mm 
focal length lenses) into a 1 mm core diameter, 20 m 
long optical fiber. The light exiting the fiber is 
collimated by lens L3 (145 mm focal length) and split 
into two paths with an uncoated optical flat (BS1). 
About 10% of the light is reflected and focused by lens 
L4 (85 mm focal length) onto PMT 1 (quantum 
efficiency - 25 %). This signal is proportional to gas 
density. The light transmitted by the beamsplitter is 
directed through a planar mirror Fabry-Perot 
interferometer (70 mm dia. mirrors, 90 % reflectivity, 
10 GHz free spectral range (FSR), finesse - 15). The 
light exiting the interferometer is focused by the fringe 
forming lens, L6. This lens consists of a pair of lenses 
(f/2 135 mm focal length and f/1.2 50 mm focal length) 
that has an effective focal length of 1600 mm. 

At the focal plane of the fringe pattern, the 
light is divided into two parts by a pair of mirrors 
(image dissector). A small mirror (6 mm dia.) is 
centered on the fringe pattern and directs light from the 
central part of the fringe through lens L7 (100 mm focal 
length) to PMT 2. This small mirror is mounted on a 
larger mirror (25 mm dia.) that directs light from the 
outer part of the fringe through L8 (100 mm focal 
length) to PMT 3. The mirrors are tilted +3° with 
respect to the optical axis. A typical calculated image 
of the inner fringe of Rayleigh scattered light is shown 
in figure 2. Note that a flow in the direction of the K 
vector, which here is in the jet flow direction (fig. 1), 


results in a positive frequency shift and increasing 
fringe diameter, while a flow in the direction opposite 
the K vector results in a negative frequency shift and 
decreasing fringe diameter. Thus, as the frequency of 
the Rayleigh scattered light increases, less light is 
detected by PMT 2 and more light is detected by PMT 
3. At the mirror location, the diameter of the image of 
the optical fiber is 11 mm. Additional optics were 
included to provide a reference image of light at the 
unshifted laser frequency. To accomplish this, several 
components could be placed in the optical path using 
remotely controlled pneumatic actuators. When placed 
in the beam path, mirror M directed laser light onto the 
diffuser (D) that scattered light into the optical fiber. 
Also, a prism assembly (PA) could be placed in the 
light path between the Fabry-Perot interferometer and 
the fringe forming lens. This served to direct light into 
a standard video camera. The video signal from this 
camera was digitized by a PC frame grabber card. A 
computer program analyzed this image and generated 
signals to control the Fabry-Perot mirror alignment. 

Uncertainty analysis 

We can now calculate the Cramer-Rao lower 
bounds given by equation 4 for velocity measurements 
based on the optical configuration described in the 
previous section. This allows us to conduct parametric 
studies to determine the optimum configuration for the 
Fabry-Perot interferometer and for the light detection 
system. The expected number of photons detected in 
time interval At by detector 1 (PMT 1) is simply RoN R At 
where Rb is the reflectivity of the beamsplitter BS and 
N r is given by equation 3. The expected number of 
photons detected by PMT 2 and PMT 3 are 

<N Dq > = (l-R B )N R AtjjS R (f)I FP (f,d r )f L 2 dfdA ( 6 ) 


where the integrations are over frequency and the area 
of the q th detector (i.e., PMT 2 or 3), f L is the focal 
length of the fringe forming lens, and I FP is the Fabry- 
Perot instrument function given by 10 


I FP {v)~ 



(7) 


where y/ is the phase change (neglecting any phase 
change on reflection) of the light between successive 
reflections given by 


c 


( 8 ) 
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Here, fi is the refractive index of the medium in the 
Fabry-Perot cavity, d is the Fabry-Perot mirror spacing, 
9 r is the angle between the ray and the optic axis, and 
F=l/(sin 2 (7i/2A r f ) where N E is the effective finesse. In 
general, the image of a monochromatic extended source 
located in the object plane consists of a series of 
unequally spaced concentric rings. In this work, 
however, the field of view is restricted by the diameter 
of the optical fiber and includes only the inner fringe as 
shown in figure 2. 

We can readily calculate the lower bounds for 
density, velocity, and temperature measurement 
uncertainties based on the detected light using three 
PMT’s as shown in figure 1. The uncertainties are 
evaluated as a function of the zero velocity fringe order 
n 0 and the flow velocity. 

It is convenient to describe the fringe location in 
terms of the fringe order rather than fringe radius because 
of the nonlinear nature of the spectral response of the 
Fabry-Perot. For example, if there is a bright fringe on 
axis, concentric bright fringes occur at integral values of 
lower orders, but the change in fringe radius decreases 
with decreasing fringe order. 

Here we define the fractional order of the fringe 
with radius r as 



(9) 


For example, at the location of the image dissector, the 
fringe order of the 6 mm diameter mirror is 0.10. 
Likewise, the fringe order corresponding to the diameter 
of the image of the optical fiber is 0.33. And, if unshifted 
laser light generates a bright fringe with radius r 0 , we 
refer to this as order n 0 . Note that here we are using the 
term ‘Iringe order” to denote change of the fringe order 
from the actual fringe order on the optical axis, which is 
2 d/?L 

The change in fringe order with optical 
frequency is 



FSR 


so the change of fringe order An with change of velocity 
AV is given by 

An = ( 2 a)vn<ej 2 ) AV ( 11 ) 

FSR 

where 8 S is the scattering angle. For example, with a 
velocity of 100 m/sec, a FSR of 10 GHz and 90° degree 
scattering, the change in the fringe order An is 0.0266. 

The lower bound for velocity uncertainty was 
numerically evaluated as a function of the fringe order 
n 0 and the velocity using equations 1, 3, 4, and 6. For 
these calculations, either two or three unknown 
parameters were assumed (i.e., the cc* in eq. 5 were N R , 


and or N R , u h and 7). The results are shown in 
figure 3 for the optical system described above for a 
single measurement using 0.5 mJ of laser energy (i.e, a 
measurement made in a 100 (usee time interval using a 
5 W laser). Note that the minimum uncertainty for 
velocity measurements depends on the fringe order n o 
and velocity. Thus, the fringe order should be set to 
minimize the velocity uncertainty for the range of 
velocities to be measured. The mean number of 
photoelectron counts in a 100 jasec time interval for the 
three PMTs are shown in figure 3a. Two cases of the 
lower bound of velocity uncertainty are shown for 
temperatures ranging from 150 to 300 K, where the 
density is fixed. One is for a flow with unknown 
temperature (fig. 3b), and one is for a flow with known 
temperature (fig. 3c). Note that, as expected, the 
measurement uncertainty is lower for colder gas. Also, 
in general, knowing the temperature significantly 
decreases the predicted measurement uncertainty. In 
this example, the uncertainty for a flow with known 
temperature is about 25 m/sec for a 300 K gas. As 
discussed above, an ideal measurement taken for a 
100 jisec time interval would have an uncertainty of 
about 4 m/sec. Thus, for our FPI based measurement 
system, the expected velocity uncertainty is about a 
factor of six larger than the uncertainty that would be 
obtained with an ideal instrument. Finally, an important 
observation is that a low velocity measurement 
uncertainty can be expected over a 600 m/sec velocity 
range (for the more general case of unknown 
temperature). 

EXPERIMENT 

Setup 

The optical system described above and shown in 
figure 1 (with the laser output power equal 5 W) was 
used to measure the flow in a subsonic free jet (8.8 mm 
exit diameter). The probe volume was located in the 
potential core of the jet to avoid any dust particles in the 
ambient air. The laboratory compressed air supply used 
for the jet was passed through a three-stage filter to 
remove particles and oil vapor. As shown in figure 1, 
the jet axis bisected the incident and scattered light 
directions. With this configuration, the measured 
velocity component was the jet axial component. Since 
the flow velocity was in the same direction as the K 
vector, the frequency shift was positive. Thus the 
Fabry-Perot interferometer fringe diameter of the 
Rayleigh scattered light increased as the flow velocity 
increased. An electronic pressure gauge was used to 
measure the total pressure in the nozzle plenum. These 
pressure measurements were used to calculate the flow 
velocity and temperature using the isentropic flow 
relations. A high precision electronic thermometer was 
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used to monitor the temperature within the FP1 
enclosure to correlate interferometer stability with 
temperature. 

Data Acquisition 

A flow chart and block diagram of the PC 
based control and data acquisition system developed for 
this work is shown in figure 4. 

The photoelectron pulses from the three PMT’s 
were amplified (Gain = 5) and sent to constant fraction 
discriminators (CFD). The CFD’s output 10 ns wide 
TTL level pulses that were counted by a PC counter- 
timer board. Pulses could be simultaneously counted on 
the three channels at rates to 80 MHz for indefinite 
times. Typical photoelectron count rates for this work 
were on the order of 5 MHz. The accumulated counts 
on each channel were recorded at preset time intervals, 
typically 100 jusec, and stored on the computer hard 
disk. The number of counts in each time interval was 
given by the difference between adjacent values of the 
accumulated counts. 

The flow chart of the software used to control all 
aspects of the instrument system operation is shown in 
figure 4a and the block diagram of the hardware is 
shown in figure 4b. The system had two modes of 
operation: (a) the Reference mode, where the automatic 
Fabry-Perot interferometer stabilization is active; and 
(b) the Rayleigh mode, where the Rayleigh scattering 
data is acquired. When the program is run, the 
sequence of operations shown in figure 4a is executed. 
The operator has the option of either switching to the 
Reference mode between each data set or taking 
sequential data sets without switching to the reference 
mode. The advantage of using the Reference mode 
between data sets is that any drift in the interferometer 
is minimized. 

Data processing 

The gas density can be obtained in a 
straightforward manner, since it is proportional to the 
total Rayleigh scattered light measured by PMT 1. 
Calibration constants (slope and intercept) are found 
from a linear least squares fit (shown on fig. 5) of the 
PMT 1 counts plotted as a function of the gas density. 

One approach to processing the PMT signals 
in order to obtain the velocity would be to use 
parameter estimation techniques based on the model 
developed for the calculation of the lower bounds of 
measurement uncertainty. This would entail doing a 
nonlinear fit of the data to the model function given by 
equation 6. However, because of the large quantity of 
data (acquired at a 10 kHz rate for the three PMT’s) a 
simpler procedure was used. This procedure took 
advantage of the fact that the temperature only varied 
over a small range, from about 295 K for low flow to 
about 255 K for the maximum flow velocity of about 
280 m/sec. The operating region was selected so the 


count rate for the inner and outer regions were 
approximately linear with respect to velocity. The first 
step was to obtain a set of mean data for the velocity 
range of interest, where the velocity was calculated 
from the total temperature and pressure ratio using 
isentropic flow relations. A linear fit was then done to 
obtain the four constants Al, Bl, A2, and B2 as shown 
for a typical data set on figure 6. 

I=A x +B x V y 0 = A 2 + B 2 V (12) 

These are solved to find the velocity as a function of the 
ratio R = O/l . 

y = A i~ A t R (13) 

B { R — 

This rational function of V = V(R) could also have been 
obtained by directly doing a nonlinear fit, but this would 
have required an iterative algorithm. 

After the calibration constants are obtained, the 
data may be processed by using equation 13 for each 
data point to obtain a velocity time history. The power 
spectrum of the velocity fluctuations can also be 
calculated. Because of the high noise levels, it is 
necessary to use relatively long data records and 
process them using a technique such as described in the 
following section. 

Welch method of modified periodograms 

We used the Welch method of modified 
periodograms 11 to calculate an estimate of the power 
spectrum of the velocity fluctuations. In this procedure, 
a long data record sampled at rate f s for time T (total 
samples = N = Tf s ) is subdivided into a number K r of 
smaller records (which may be overlapping), each of 
length L samples. The modified periodograms of each 
sub-record are calculated using a data window; these 
individual periodograms are then averaged to obtain the 
estimate of the power spectrum. The frequency 
resolution of the resulting spectrum is thus f s /L. By 
overlapping the segments by one half of their length, a 
near maximum reduction in the variance in the spectral 
estimate is achieved; the variance in the estimated 
spectrum is reduced by a factor of 11/9 K r compared to 
the variance of the spectral estimate calculated directly 
from the original long data record. 

RESULTS 

Although the eventual application of this 
technique is the study of turbulent flows, the 
measurements reported here were all taken in the 
potential core of the subsonic free jet where the 
turbulence level should be quite low. This initial work 
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in the free jet allowed the system to be evaluated for 
known flow conditions, where the flow velocity, 
density, and temperature could easily be determined 
using isentropic flow relations. 

An example of a set of mean velocity 
measurements (each calculated using 100,000 samples 
taken at a 10 kHz rate) is shown in figure 7. The solid 
line is the velocity calculated from the measured total 
and ambient pressures using isentropic flow relations. 
The calibration constants described in the data 
processing section were obtained from data taken five 
days earlier. The rms deviation of the mean velocities 
from the isentropic values is 7.7 m/sec. This plot shows 
that the Fabry-Perot stabilization system is able to 
reproduce the correct interferometer alignment (locking 
to a preset fringe diameter before each measurement). 

The distribution function of the counts from 
PMT 1, which are proportional to density, is shown for 
two flow velocities (22 and 284 m/sec) in figure 8. The 
square of the standard deviation of the distribution (the 
variance) is very close to the mean, as expected for the 
Poisson statistics of the photoelectron counts. Also 
note that at the higher velocity, the density is greater 
because the temperature is reduced by the isentropic 
expansion. 

The velocity distributions of two time history 
records are shown in figure 9 along with Gaussian fits. 
Since the measurements were taken in the potential core 
of the jet, the actual turbulence intensity is low and the 
fluctuation in the velocity is due almost entirely to the 
shot-noise. The measured values of velocity standard 
deviation, about 20-30 m/sec, agree with the predictions 
of the lower bounds calculations shown earlier. The use 
of the standard deviation as a measure of turbulence 
intensity can only be done for velocity fluctuations on 
the order of, or larger than the standard deviation due to 
the shot-noise. In this situation, since the shot-noise 
and velocity fluctuations are independent and add in 
quadrature, one can subtract the contribution of the 
shot-noise. However, we can see that this is much more 
attractive for higher velocity flows. 

Finally, two examples of velocity power spectra 
calculated using the Welch method of modified 
periodograms described above (with K = 195 

overlapping segments of length L = 1024) are shown in 
figure 10. For the low velocity flow (22 m/sec) shown 
in Figure 10a, the spectrum is relatively flat, except for 
very low frequencies. For the higher velocity flow 
(284 m/sec) shown in Figure 10b, the spectrum shows 
significant power in the frequency range 200-800 Hz. 
This is most likely caused by flow fluctuations 
generated in the nozzle air supply line. 


CONCLUDING REMARKS 
An improved technique for obtaining dynamic 
gas density velocity measurements using molecular 
Rayleigh scattering was described. An analysis was 
presented that established a lower bound for 
measurement uncertainty of about 20 m/sec for 
individual velocity measurements. The velocity was 
determined by analyzing the scattered light with a 
Fabry-Perot interferometer. The signals from three 
photomultiplier tubes were simultaneously recorded 
using photon counting at a 10 kHz sampling rate and 10 
second recording periods. Software and hardware 
interfaces were developed to allow computer control of 
all aspects of the experiment and data acquisition. 
Density and velocity measurements taken in a subsonic 
free jet were presented. The measured velocity 
standard deviation was close to the lower bound 
prediction. Velocity and density power spectra were 
obtained using the Welch method of modified 
periodograms. The stabilization system used for the 
Fabry-Perot interferometer allowed calibration data 
obtained one day to be used for data acquired five days 
later. This technique, which simultaneously measures 
time history of density and velocity, will be an 
important new tool for study of noise generation in 
turbulent flows. 
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Fig. 1 - Layout of Rayleigh scattering system 
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Fig. 2 - Fabry-Perot interferometer fringe showing two 
regions where light is directed to PMT 2 (INNER) and PMT 
3 (OUTER). 
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Fig. 3 (a) Number of counts for PMT 1 , PMT 2 
(INNER), and PMT 3 (OUTER) as function of 
velocity for zero velocity fringe order n 0 = 0.10; 
(b) Lower bound for velocity uncertainty for 
unknown temperature; (c) Lower bound for 
velocity uncertainty for known temperature. 
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Fig. 5 - Calibration plot for gas density . 
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(b) 

Fig. 4 - Experiment control and data acquisition system (a) Flow 
chart of system operation; (b) Block diagram 
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Fig. 6 - Linear fit of INNER (a) and OUTER (b) counts 
as function of velocity 
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Fig. 7 - Velocity obtained from Rayleigh 
scattering measurements and from isentropic flow 
relations. 
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(b) 

Fig. 8 - Distribution of PMT 1 counts (representing 
density) for mean flow velocity of 22 m/sec (a) and 
284 m/sec (b). 
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Fig. 9 - Velocity distribution function of mean 
flow 22 m/sec (a) and 284 m/sec (b). 




Fig. 10 - Power spectrum of velocity 
fluctuations for mean velocity = 22 m/sec (a) 
and 284 m/sec (b). 
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